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^: ABSTRACT 

i-^ ' In the last few years considerable progress has been achieved in our understanding of the 

decays of heavy flavour hadrons. One can now calculate inclusive transition rates in QCD 
k>( \ proper through an expansion in inverse powers of the heavy flavour quark mass without 

;_( ■ recourse to phenomenological assumptions. The non-perturbative contributions are treated 



C^ 



systematically in this way; they are found to produce corrections of order a few percent in 
beauty decays, i.e. typically somewhat smaller than the perturbative corrections. One finds, 
among other things: (a) The lifetime of B~ mesons is predicted to be longer than that of B^ 
mesons by several percent, (b) The QCD prediction for the semileptonic branching ratio of 
B mesons appears to exceed present experimental values. We discuss the implications of this 
discrepancy. The phenomenological engineering that has been developed for the description 
of exclusive two-body modes of B mesons has reached a mature stage and awaits more 
precise and detailed experimental tests. First steps towards a genuine QCD treatment of 
these modes are being made. 



'Anyone who keeps the ability to see beauty never grows old' Franz Kafka 



^To appear in the second edition of the book 'B Decays', S.Stone (ed.), World Scientific 
^Permanent address 



1 Introduction 

In the two years since the appearance of the first edition of this book our under- 
standing of the decays of beauty hadrons has been improved quite significantly. The 
availabihty of better and more comprehensive data has obviously helped here - but 
so has the emergence of more reliable theoretical tools! Rather than having to rely on 
phenomenological models we can now employ approaches that are directly based on 
QCD with no appeal to a 'deus ex machina', at least for semileptonic and inclusive 
non-leptonic decays. Therefore the article on non-leptonic beauty decays had to be 
re-written completely. Inclusive decays of beauty hadrons can now be treated in QCD 
in a quantitative fashion with systematic estimates of the uncertainties. 

The employed formalism combines the heavy quark mass expansion with other 
elements derived from QCD proper. Its technical foundation is the Wilson operator 
product expansion (OPE) [|l|. The idea of applying OPE to the inclusive heavy flavour 
decays had emerged in in the eighties and has now grown into a well-developed 
scheme §1,1,11. 

The basic procedure can be illustrated by a simple analogy with nuclear jS decay. 
There are two effects distinguishing leptons in the decay of a neutron bound inside a 
nucleus from those in a decay of a free neutron: 

(a) nuclear binding effects; 

(b) Pauli statistics correlating the electrons surrounding the nucleus with those 
emerging from the (3 decay. 

The typical energies of the bound electrons e^i are certainly small compared to 
Erei, the energy released in the decay; let us assume - although this is not true in 
reality - that also the nuclear binding energies enud were small compared to Ej-ei- In 
that case nuclear (3 decays would proceed like the decays of free neutrons to a good 
approximation; corrections to this simple 'spectator' picture could be computed via 
an expansion in powers of enud/Erei, ^ei/Erei In practice, however, the corrections for 
nuclear (3 decay are incorporated by explicitly using the wavefunctions of the bound 
nucleons and electrons. 

There arise analogous corrections to the decay rate for a quark Q inside a hadron: 

(a) interactions of the decaying quark with other partons in the hadronQ. 

(b) Pauli interference effects of the decay products with other partons in the 
hadron; e.g.: bu — >• cudu. 

The difference with the example above is quite obvious: even in the limit mq —>■ oo 
a non-relativistic bound state treatment is inapplicable since the dynamical degrees 
of freedom of the heavy flavour hadron Hq cannot be fully described by a hadronic 
wavefunction! We will return to this point later on. The most reliable approach is 
then to evaluate weak decay rates of heavy flavour hadrons through an expansion in 
powers of fihad/^TT'Q] "^Q ^ the heavy flavour quark mass - is a measure of the energy 
release in the decay and fihad represents ordinary hadronic scales which enter through 
the matrix elements for heavy flavour hadrons. Since fihad < 1 GeV (details will be 
given below) one expects such an expansion to work quite well for beauty decays. 

^There is also the annihilation of the heavy quark with the light (spectator) antiquark for which an 
analogy is found in the K capture of bound electrons by a heavy nucleus. 



The example given above illustrates two important features of our analysis: it 
applies to inclusive non-leptonic and semileptonic transitions, and the usefulness of 
such an expansion rests on a large energy release in the decay; i.e. in Qi — >■ Q2W\ 
where W denotes q'q or lu one requires mq^ — mq^ ^ fJ^had- This is in contrast to 
the approach usually referred to as the heavy quark symmetries which applies to 
exclusive semileptonic amplitudes provided that the initial and the final quarks are 
both heavy, i.e. mq-^, mq^ ^ ^^ad but the ratio mq^/mq^ is otherwise arbitrary. The 
QCD-based description of the inclusive decays of heavy flavours naturally incorporates 
many elements of Heavy Quark Effective Theory (HQET) [§, ^. The systematic 
expansion in rriq^ is inherent to both, and it allows us - as we will see shortly - to 
determine the size of important input parameters for our analysis. 

This survey will be organized as follows: in Sect. 2 we review the general features 
of the l/mq expansion for inclusive decay rates; in Sect. 3 we present quantitative 
predictions for lifetime ratios, semileptonic branching ratios and radiative widths and 
discuss their theoretical uncertainties; in Sect. 4 we address exclusive non-leptonic 
decays before giving a summary and an outlook in Sect. 5. 

2 Inclusive Decay Rates of Beauty Hadrons: Formalism 

As mentioned above, the systematic approaches to inclusive heavy flavor decays date 
back to the beginning of the eighties. Unitarity relates inclusive decay rates to the 
imaginary part of certain forward 'scattering' amplitudes; this is a rather trivial 
observation ||l^, yet it opened the way for the consistent use of OPE. The starting 



object of our analysis is the transition operator T{b ^ f —>■ b) describing the forward 
scattering of b quarks via an intermediate state /. To second order in the weak 
interactions the transition operator is given by [|I^ 

f{b^f^b)=ij d^x{Cw{x)C^y,mT. (1) 

where Cy/ denotes the relevant effective weak Lagrangian and {.}r is the time-ordered 
product. Treating mi, as a large parameter a Wilson OPE allows expressing the non- 
local operator T as an infinite sum of local operators of increasing dimension with 
coefficients containing higher and higher powers of 1/mf,. The lowest dimensional 
term will dominate in the limit nih -^ oo; for beauty decays that is the dimension 
three operator bb. The width for the decay of a beauty hadron Hi, into an inclusive 
final state / is obtained by taking the expectation value of T between the state Hh. 
Through order 1/mf one finds [|, |, |], 0: 






{H,\bb\H,) cif) {H,\Ua^,G,,b\H,) 






(2) 



where the dimensionless coefficients Cj(/) depend on the parton level characteristics 
of / (such as the ratios of the final state quark masses to rrib), KM denotes the 



appropriate combination of weak mixing angles and G^i, the gluonic field strength 
tensor. The last term implies also the summation over the four-fermion operators 
with different light flavours q. Notice that the factor 1/2Mh^ reflects the relativistic 
normalization of the state \Hb). 

It is through the expectation values of local operators appearing on the right-hand 
side of eq. (2) that the dependence on the decaying hadron, and on non-perturbative 
forces in general, enters instead of entering through the wavefunctions as in nuclear f3 
decay. Since these are matrix elements for real b hadrons one sees that T{Hi, -^ f) is 
indeed expanded into a power series in Hhad/f^h- The heavy quark expansion enables 
us to express decay rates in terms of the expectation values of local operators taken 
between beauty hadrons. Using heavy quark expansions one can relate some of these 
hadronic expectation values - and in particular those that appear in the leading terms 
of the expansion in eq. (2) - to other observables and thus extract their size, as we 
will discuss now. 

(i) To determine {Hi,\bb\Hi,) one expresses, via equations of motion, the operator bb 
through another series in inverse powers of mi, (which constitutes a static expan- 
sion) P: 

{H,\bb\H,) = 

mv,h,b - -^mv ■ Df - {iDf]b + ^bia ■ Gb\H,) + 0{l/ml) (3) 

with v^ denoting the four- velocity of the heavy hadron. The first operator appearing 
on the RHS of eq. (3) is actually the Noether current for the (global) heavy flavour 
quantum number; its expectation value is thus determined by the beauty content of 
Hh and therefore 

{H,\bb\H,)/{2MH,) = 1 + 0{l/ml). (4) 

It is this term that exactly reproduces the parton model spectator result in the limit 
rrih -^ oo, which attributes equal lifetimes to all hadrons of a given heavy flavour, 
(ii) The chromomagnetic operator bia ■ Gb appears directly in eq. (2) and indirectly 
through the expansion of bb. Its expectation value vanishes for the baryon A;,: 

{Ab\bia ■ Gb\Ab) - (5) 

For the B meson it is given by the hyperflne splitting of the B* and B masses: 

{B\bi(T ■ Gb\B)/{2MH,) ^ -(M|. - M|) ^ 0.74 GeV ^ (6) 

(iii) The second operator on the RHS of eq. (3) describes the kinetic energy of 
the b quark moving under the influence of the non-trivial gluon background fleld 
prevalent inside the hadron Hb: {Hb\b[{iv ■ Df - {iDf]b\Hb) - {Hb\b{iDfb\Hb) = 
{{pby)Ht ■ (2Mi:/J. Its appearance in eq.(3) has a very transparent meaning: 



2Mh. 2ml 8 m 



1 - ' 2 + r 2 + 0{l/ml) (7) 



'H, ^"ife 



The flrst two terms on the RHS of eq. (7) represent the mean value of the factor 
Vl — 'ip' due to time dilation slowing down the decay of the b quark in a moving 



frame. The heavy quark expansion relates the difference in the expectation values of 
the kinetic energy operator for heavy flavour baryons and mesons to the masses of 
the charm and beauty baryons and mesons |12[ : 



mr)A, - mr)B ^ J^^ ■ {ms) - m^j - [(m^) - m^j} (s) 

where {Mb^d) denote the 'spin averaged' meson masses, 

{Mb) = \{Mb + ?>Mb'^) (9) 

and likewise for (Mo)- Eq. (9) implies that c quark can be also treated as heavy so 
that {{pb)'^)Ht = {{Pc)'^)hc- From the present data we obtain 

{iPb?)A, - {{Pbf)B - (-0.03 ± 0.17) GeV ' (10) 

where the error is due to the ±30 MeV experimental uncertainty in M\^. The A^ 
mass has to be measured with better than a 10 MeV precision to make this relation 
numerically useful. 

Using the commutator algebra of the covariant derivatives iD^^ one can derive an 
'uncertainty principle' for their components and thus establish a model independent 
lower bound p3[| : 

mf)B> 0.18 {GeV f. (11) 

An existing analysis based on QCD sum rules yields a value only three times larger 
than this lower limit [Q: 

{{p,f)B^ 0.6 {GeV f. (12) 

(iv) The expectation value for the four-quark operators looks very similar to the one 
controlling B^ — B^ oscillations: 

{B{p)\{hLl,qL){qLlM)\B{p)) ~ flp^p. (13) 

where we have set the so-called bag factor to unity. 

Before discussing the phenomenology that is obtained from eq.(2) we want to point 
out seven basic observations: 

(1) The most important aspect - quantitatively as well as conceptually - of the 
expression in eq. (2) is contained in the element that is missing there: there are 
no non-perturbative contributions of order l/rrift to fully integrated rates [^! The 
numerical impact of this fact is obvious: since the leading non-perturbative corrections 
then arise on the 1/m^ level, they fade away quickly with increasing heavy flavour 
quark mass and amount, for the beauty decays, to effects of order of several percent 
only since the scale is set by the quantities 

{B\bia-Gh\B) , , 

^-- 24-(2Mi) -°-°^^' ^^^«) 

Kb = -^^ ~ 0.015. (146) 

mi 



Let us note in passing that the analogous contributions in charm decays are much 
larger since they are amplified by a factor {mb/rric)'^ ~ 10! 

The conceptual relevance of the absence of l/mt, terms is of a more subtle, but not 
less important nature. On the one hand it confirms the conjecture that the (current) 
quark mass rrih rather than the hadron mass Mb represents the natural expansion 
parameter. For if the total width were correctly expressed in terms of Mb then a term 
linear in l/rrih had to appear: T{B) oc G'^M^ = G'jpml{l + 5As/m(, + ...) with the 
notation Mb = rrih + A^ + 0{l/mf,). This is actually more than an academic point 
for it would have a significant impact on the lifetime difference between B mesons 
and Afe baryons: (r(Afe) — t{B))/t{B) oc l/rrib since Aa^^ 7^ A^! On the other hand 
it is quite instructive to understand the dynamical reason why no linear l/m?, terms 
arise in total decay rates whereas they do appear in mass formulae. Looking at the 
explicit QCD calculation given in ref. one can already infer why no contribution 
of order l/irti, arises: only a dimension four operator could generate such a term and 
there simply does not exist such an appropriate operator that is gauge invariant and 
cannot be absorbed into the dimension three operator hh by means of equations of 
motion. (This means that proper care has to be applied in employing a consistent 
definition of the heavy fiavour mass m;,.) It is therefore the colour symmetry, i.e. the 
fact that the colour fiow is conserved that ensures the absence of a l/r/ife correction. 
This correction is absent also for differential distributions like the lepton spectra in 
semileptonic decays outside the end-point region. The size of this end-point domain 
is of order f^hadr/inb if all energies are measured in units of m^. In this domain we 
cannot limit ourselves to the operators of the lowest dimension, and one needs to 
sum up an infinite series to determine the shape of the distribution which leads to 
modification of the spectrum of order unity in the end-point region. It is remarkable 
that the full integral over the end-point region still has no corrections proportional 
to l/rribl 

(2) The distinction between beauty baryon and meson decays is systematically 
expressed through differences in the appropriate expectation values of the same op- 
erators, see eq.(2). Differences arise first in the leading non-perturbative corrections 
of order 1/m^, see eqs.(5, 6, 8). Yet apart from some small SU{3)fl breaking effects 
they affect B^, B^ and Bg decays in the same way; likewise for A^ vs. D^, D^ and 
Dg decays. Numerically they are comparable to - actually typically smaller than - 
perturbative corrections in beauty decays whereas they dominate perturbative effects 
in charm decays. 

(3) Local four-quark operators of dimension six finally produce differences between 
all the B meson lifetimes: 

^ nonspecty^ ) JB 

Tspect{B) rrif^ 

which formally scales like 1/m^. These effects are therefore predicted to be greatly 
reduced relative to the considerable lifetime differences in the D system. 

(4) Contributions of order l/mf are generated by dimension seven operators. Yet 
it appears to be practically unfeasible to determine the expectation values for all or 
even most of them. What seems possible - although it has not been done yet - is to 
analyze a small subset of them, namely those yielding factorizable contributions, for 



a more detailed error estimate. 

(5) It is intuitively obvious that the b quark does not rest inside a beauty hadron 
Hb, but will move around with a certain 'Fermi momentum'. This has been incorpo- 
rated into phenomenological models of inclusive heavy flavour decays, first in ref. |T^ . 



It has been first stated in ref. |0 and then further discussed in refs. [0, |T^ |53|, ^ 
that this notion of Fermi motion of a heavy flavour quark finds a natural home also 
in a rigorous QCD treatment; yet, strictly speaking, no hadronic wavefunction for Hh 
in the usual sense can be found that reproduces the Fermi motion beyond the second 
moment {{p}^)- We will not discuss this in any detail here since {{p)"^) has a reduced 
numerical relevance in fully integrated rates (although it is of crucial importance for 
shaping the end-point spectrum in the inclusive semileptonic and radiative decays). 

(6) While our expressions for inclusive decay rates are firmly based on QCD, one 
has to invoke explicitly - and at the moment additionally - the concept of duality. 
The operator product expansion is unambigously defined in the Euclidean domain, yet 
the kinematics of the actual decay are of a time-like Minkowskian nature. It is then 
conceivable - although it has never been illustrated in a clear way - that a translation 
between the Euclidean and the Minkowskian expression that is based on local duality 
does not hold in non-leptonic or for that matter even in semileptonic decays. We view 
this as a mathematical rather than as a physical caveat. The conjecture of duality 
can be supported by some general arguments, yet their discussion would lead beyond 
the scope of this article. 

(7) The expression in eq. (2) is based on OPE where one separates short distance 
and long distance dynamics by isolating the latter in the local operators and their 
matrix elements and letting the former determine the c number coefficients Cj(/). In 
actual calculations one goes one step further: one computes the coefficients Cj(/) in 
perturbation theory alone although non-perturbative short distance contributions do 
exist. The latter are guestimated (and in r decays found ||TB[) to be quite small. We 



adopt this procedure which we refer to as the 'Standard Version' of OPE although 
we will also comment on possible limitations later on. 

3 Phenomenology of Inclusive Beauty Decays 

There are five types of inclusive observables we will discuss here, namely: (1) total 
lifetimes; (2) semileptonic branching ratios; (3) other inclusive non-leptonic decays; 
(4) radiative decays; (5) charm multiplicity in the final state and (6) B -^ charmonia 

+ X. 

3.1 Total Lifetimes 

There exists a triple motivation behind measuring the lifetimes of different species of 
beauty hadrons as precisely as possible: (a) Representing the most inclusive quantity 
it provides a clear and well-defined challenge to theory, (b) It allows to obtain the 
semileptonic width from the measured semileptonic branching ratios; from this width 
one extracts the KM parameter \Vcb\ etc. (c) It is a pre-requisite of a detailed analysis 
of B^ — B^ oscillations. 



There is no basic uncertainty about the weak forces driving non-leptonic beauty 
decays: at the scale My/ they are given by the Lagrangian 

£^^=i(/i = Mw) = -^[VcCLi^hL + K^^lTm&l] ■ K^dLl^^UL + V:-SLlM (15) 

where we have ignored Cabibbo suppressed transitions and also the b -^ t coupling 
since we will not discuss B^ — B^ oscillations and Penguin contributions here. Ra- 
diative QCD corrections lead to a well-known renormalization at scale m;,, which is 
often referred to as ultra-violet (UV) renormalization: 

^W^^^if^ = ^b) = —7=-VchV*d{Cl{cLl^J:)L){dLli,UL) + C2{dL'^^,hL){cL'l fiUL)} (16) 

for h -^ cud and likewise for h -^ ccs etc. transitions; the QCD corrections are lumped 
together into the coefficients ci and C2 with 

ci = 2(c+ + c-), C2 = -(c+-c_) (17a) 



as[m^) 33 — 2Nf 2 

in the leading log approximation. Numerically this amounts to 

Ci(LL) ~ 1.1, C2(LL) ~ -0.23 (18) 

Including next-to-leading log corrections one obtains 

ci(LL + A^LL) ~ 1.13, C2(LL + A^LL) ~ -0.29 (19) 

i.e. a mild enhancement of the original coupling together with the appearence of 
an induced operator with a different color flow. Later we will also include [1^ the 



so-called 'hybrid' renormalization reflecting radiative corrections in the domain from 

TTlb down to fihad [H • 

As already stated in Sect. 2 differences between B meson lifetimes arise on the 
1/ml level generated by local four-quark operators (bL'lfi(lL){(lLlubL)- Based on this 
scaling law one can already infer from the observed D meson lifetime ratios that the 
various B meson lifetimes will differ by no more than 10 percent or so. 

In phenomenological models two distinct mechanisms producing lifetime differ- 
ences had been noted, namely 

• Weak Annihilation (WA) and 

• Pauli Interference (PI) |^ 

in qualitative analogy to the situation in nuclear jS decay as explained in the In- 
troduction. In the l/rrifj expansion they emerge as follows. There are two types of 
four-quark operators which are distinguished by how the light quark flavours are con- 
nected inside the hadron Hb. This can be seen from Figs. 1. Upon integrating out the 
c, q and q' fields in the diagram of Fig. la where the square boxes represent C^P^^ 
one obtains the operator bb; cutting then the q' line in Fig. la and connecting it to 



the q' constituent of the B meson, as shown in Fig. lb, one has a WA transition 
operator; cutting instead the q hne and connecting it to the B constituents, see Fig. 
Ic, leads to the four-fermion operator describing PI. 

It turns out that the WA processes can change B lifetimes by no more than, say, 
1%; due to interference with the spectator reaction they could actually prolong T{Bd) 
relative to t{B~) rather than reduce it ^! The dominant effect is provided by PI 
which produces an additional contribution to the B~ width: 

T{B-) = Ts,ect{B) + ^Tpj{B-) (20a) 

AFpKi?-) - To ■ 24vr^^[4 - cl + ^(4 + ^)], T, = ^\V{cb)\^ (206) 

where the appearance of fs reflects the fact that PI - like WA - requires the spatial 
overlap of two (anti) quark fields. 

Eq. (20b) exhibits an intriguing result: AFp/(i?~) is positive for c+ = 1 = c_, i.e. 
PI acts constructively. Radiative QCD corrections with c+ ^ 0.84, c_ ~ 1.42 turn 
PI into a destructive interference which prolongs t{B^) by a tiny amount. In eq. 
(20b) only UV renormalization has been applied. Hybrid renormalization amplifies 



this effect considerably and one obtains H, 19 



f 2 r'^ 4- r'^ 1 



ATpjiB-) ^ Fo ■ 2An^j^n-%ci - ci^^ + ^-^^ - ^-{n^r^ - 1)(4 - c^_)]. 



K^P^Yl\h=n-\n, (21) 



Altogether one finds: 



f2 



^^^-■^'0.05-- {f-^TT^, (22) 



T{Bd) (200 MeV 

i.e. the lifetime of a charged B meson is predicted to exceed that of a neutral B 
meson by several percent, but not more than ten percent [[ Corrections of order 1/m^ 
which have been ignored here are unlikely to change this prediction significantly; this 
statement will however be qualified below in our discussion of the semileptonic B 
branching. 

One also expects 

f{Bd) - f{Bs) (23) 

where f denotes the average lifetime of the two mass eigenstates in the B^ — B^ 
system. It is at least amusing to note that the largest lifetime difference among B 

^It should be noted that an analogous analysis yields t{D^)/t{D'^) ~ 2, which is quite consistent with 
the observed value of ~ 2.5. Yet one has to keep in mind that the charm quark mass is not much larger 
than typical hadronic masses; the convergence of the l/rric expansion is thus too slow, if it happens at all, 
to yield a better than semi-quantitative description of non-leptonic charm decays. For a recent review of the 
theoretical situation here see ref. ESfl. 



mesons is most likely produced by a subtle mechanism, namely Bg-Bg oscillations 
with both Pl-like and WA-like mechanisms contributing [p3|: 



AT{Bs) _ T{Bs,short) - T{Bs^iong) _ ^ IS J'' 



r{B,) ~ r{B,) ~ °'^^ ■ (200 MeV )2 " ^^^^ 

One can search for the existence of two different Bg lifetimes by comparing t{Bs) as 
measured in Bg -^ ipcj) and in Bg -^ IvX. Analogously one can compare t{Bci) as 
obtained from B^ — > ^-f^s, from Bfi -^ ipK* and from B^ — > IvX. In the B^ case 
one theoretically expects a lifetime difference on the percent level only. Whether an 
effect of the size indicated in eq. (24) is large enough to be ever observed in a real 
experiment is of course a different matter. It has to be said, though, that eq.(24) 
does not represent a 'gold-plated' prediction. It is conceivable that the underlying 
computation underestimates the actual lifetime difference. 

No detailed analysis has been performed yet on t(A6); simple estimates lead to 
the expectation 



r(A 



b) 



<Ba) 



0.9. (25) 



Present measurements yield 

r(fi+) = 1.66 ± 0.11 psec (26a) 

r(5rf) = 1.51 ±0.10 psec (266) 

\^! = 1.12 ±0.09 (26c) 

T{Bg) = 1.54 ± 0.24 psec {26d) 

T{Ab) = 1.07 ± 0.16 psec (26d) 

While the predictions stated above on the lifetime ratios for B mesons, see eqs. (22,23), 
are quite consistent with the measurements, one cannot draw a definite conclusion at 
the moment. The At, lifetime appears to be shorter than the B lifetimes, though: 

r(fi+) -r(Ab) = 0.59 ±0.19 psec (27) 

r(Srf)-r(Af,) = 0.44 ±0.19 psec (28) 

Qualitatively this is expected although the size of the difference seems to be larger 
than anticipated. 

3.2 Semileptonic Decays 

The semileptonic branching ratio of beauty quarks depends on fundamental quanti- 
ties of the Standard Model, namely KM parameters and quark masses. It is then 
the primary goal of the measurements to determine the numerical values of these 
quantities from the data. The theoretical challenge on the other hand consists in 
disentangling the effects of the strong interactions both in their perturbative as well 
as non-perturbative aspects. 

10 



The recent experimental studies have reached a new level of accuracy and relia- 
bility: a 'model-independent' ARGUS analysis yields [ pS] 

BRsl{B) = 9.6 ± 0.5 ±0A% (29a) 

whereas the CLEO collaboration finds ^^ 

BRsl{B) = 10.65 ± 0.05 ± 0.33% (296) 

using the model of Altarelli et al. for the shape of the lepton spectrum. It is an 
intrinsically phenomenological description, yet one should keep in mind that it pro- 
vides a practically good approximation to the true QCD lepton spectrum as calculated 
through a l/mf, expansion 0. The present data thus clearly suggest: 

BRsl{B)U.p< 11% (30) 

In a naive parton model where even perturbative QCD is ignored one obtains 

BRsLib -^ clu) ^ 15 ^ 16% (31) 

i.e. a non-leptonic enhancement of ~ 50% has to be found to reproduce the data. At 
first sight this would not seem to represent a stiff challenge - yet so far we have failed 
to meet it! 

There are non-perturbative as well as perturbative corrections to the semileptonic 
beauty branching ratio; our ignorance about the former was very considerable before 
the arrival of the Heavy Quark Expansion; we will discuss them in sequence. 

3.2.1 Non-perturbative Corrections to BRsl 

The semileptonic and non-leptonic widths through order l/ml are given by || 



Tsl{B) = To ■ 



{B\bh\B) 



2M, 



B 



/o(x,0,0) + ^(x-^-2)Jo(x,0,0) 
mi ax 



% 



(32) 



(^l^^l^) L i^r^^. , ^^h,d 



T^B) = T,-Nc- ' ^^^ ' ■ |Ao[S/o(x) + ^(x- - 2)S/o(a;)]- 

8 ^2 ^ ■ [/2 (x, 0, 0) + /2 (x, x, 0)] I . (33) 

where the following notations have been used: /q and I2 are phase-space factors: 
Jo(x,0,0) = (l-x2)(l-8x+x2)-12xMnx, /2(x,0,0) = (l-x)^ x = {mjnikf (34) 

1 + v 



Jo(x, X, 0) = v{l - 14x - 2x^ - 12x^) + 24x^(1 - x^) In , v = Vl - 4x (35a) 

1 — V 

hix, X, 0)=v(l + - + 3x^) - 3x(l - 2x2) In ^-^ (356) 

2 1 — V 
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with /o,2(a^, X, 0) describing the h -^ ccs transition, and SJo(x) = /o(x, 0, 0)+/o(x, x, 0); 
Aq = TjJ, A2 = (c^ — c^)/6 where r] = (c^ + 2c^)/3, and J represents the effect of 
the subleading logarithms [^ (unknown for A2). Moreover, 

f^l^^{B\hza-Gb\B). (36) 

The matrix element {B\bb\B) enters as an overall factor into both the semileptonic 
and non-leptonic width; its value does therefore not affect the semileptonic branching 
ratio. On the other hand {B\bi(T-Gb\B) which is determined from the observed B* — B 
mass splitting does, and it actually reduces BRsl{B) since A2 < 0! Using 

mf'' = 4.8 GeV , mf ^ ~ 1.35 GeV (37) 

one finds 

5BRsLiB)\,,on-pert ~ "0.02 • BRsl{B) ^ -0.003 , (38) 

i.e. a very small reduction! |^ Corrections of order 1/ml have been analyzed in ref. ||27| |; 
as expected, they are estimated to be quite insignificant. This means that non- 
perturbative corrections can to a good first approximation be ignored in BRsl{B) ! 

3.2.2 Perturbative Corrections to BRsl 

The preceding discussion shows that it is mainly the perturbative corrections that 
control the size of BRsl{B). They indeed generate a sizeable non-leptonic enhance- 
ment thus reducing BRsl{B), as desired; yet numerically they fall short of the goal. 



For one finds |27 



BRsl{B)\qcd>12.^% (39) 

It has been known for some time that the measured semileptonic branching ratios fall 
below the predicted ones. Now, however, one has reached a stage where one has to 
take such a deficit seriously since both the experimental and the theoretical analyses 
have become rather mature. 

3.2.3 Fabula Docet? 

An intriguing problem has arisen, which warrants serious consideration: how can one 
find an additional non-leptonic enhancement of at least 15 to 20 % to satisfy the 
bound of eq.(30)? 

There are various scenarios for resolving this apparent puzzle: 
(i) Improved data could move BRsl{B) above 12 % . 

(ii) There could be a 'cocktail', i.e. a combination of several smallish effects all 
working in the same direction: the experimental number could inch up; higher order 
perturbative and non-perturbative corrections could turn out to be somewhat larger 
than estimated by us. 

■^It should be noted that the corresponding effect is much larger in D decays - it leads to a roughly 50 % 
reduction in BRsl{D) and is essential for a self-consistent understanding of charm decays. 
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(iii) Non-perturbative corrections could conceivably be dramatically larger than an- 
ticipated by us. This certainly would require going beyond the 'Standard Version' 
of OPE; for it would presumably mean that even in beauty decays there are numer- 
ically significant non-perturbative corrections that enter through the coefficients in 
the OPE. In that scenario one would probably obtain considerably larger differences 
in the lifetimes of beauty hadrons than stated in eqs.(22, 25) above! 
(iv) The most intriguing possibility would be the intervention of New Physics in B 
decays. 

3.2.4 Lepton Spectra 

The method outlined above has been extended to treat the lepton energy spectra 
in Hb -^ IvX transitions. The expansion is now in 1/(1 — y)mi, rather than in 
l/rub, where y = Ei/ E"^"-^ denotes the normalized lepton energy (for h -^ u decays 
y = 2Ei/'mfy). This expansion is obviously and necessarily singular at y = 1, i.e. in 
the endpoint region, and one has to apply care in interpreting the results there. 

To order 1/ml the spectrum dV/dy is evaluated without any free parameters with 
the non-perturbative corrections entering through the expectation values Gb and Kb] 
in practice there is at present some numerical uncertainty in the size of the kinetic 
energy term Kb, as already mentioned. The shape of the spectrum thus derived from 
QCD turns out to be remarkably similar to the one obtained from the phenomeno logi- 
cal AC'^M'^ model [0 that has been fitted to the data. An important element of that 
model was the introduction of a Fermi motion ascribed to the heavy flavour quark. 
It has been found [^, |1^, [1^ that this notion of Fermi motion finds a natural home 
also in a rigorous QCD treatment: one can define a universal distribution function 
that describes the motion of the beauty quark inside Hi, irrespective of the specifics 
of the decay process; to obtain the observable spectrum predicted from QCD one 
has to fold this distribution function with the primary lepton spectrum from beauty 
decays Q. The properties of this distribution function control the shape of the lepton 
spectrum in the endpoint region. The second moment of the distribution function is 
given by (p^). One can also show that this function as a matter of principle cannot be 
obtained from a non-relativistic hadronic wavefunction, although the b quark Fermi 



motion is of a non-relativistic nature 13 



Lastly, one finds that sizeable differences can arise in the endpoint spectrum of 
Bfi vs. B~ mesons to order 1/ml due to WA in the KM suppressed decays [0. A 



detailed study of this difference in the spectra can provide information about the 
four-fermion matrix elements driving WA in nonleptonic decays and affecting PI. 

3.3 Other Inclusive Non-Leptonic Decays 

The l/rrih expansion allows not only to calculate the overall non-leptonic and semilep- 
tonic rates, but also various sub-classes, like, e.g., those non-leptonic transitions that 
are driven by 6 ^ ccs, b -^ uud etc. While we have expressed some notes of caution 

^One should remember however that the universaUty holds only as long as the mass of the final state 
quark is the same. In QCD the corresponding functions that enter here are quite different for 6 — > c and for 
b —> u (or b —> s + j) transitions. 
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about the prediction for b -^ ccs reactions since the energy release there is not very 
large, no such caveat applies to the b — > uud process: the expression in eq. (33) 
is easily adapted by the obvious substitutions: Vub for Vbc and m„ = for mc, the 
expectation values of the local operators remain the same. Whether the prediction 
for the KM suppressed, inclusive non-leptonic B decays can ever be tested with any 
decent accuracy, is of course a different matter. 

3.4 Radiative Beauty Decays 

The non-perturbative contributions to T{Hh —>■ s+7) through order 1/m^ are obtained 
in a straightforward manner [EUI: 



\ 2MHi 2MH,mf J 

where we have set rus = 0. From eqs.(40) and (32) one reads off that the ratio 
r{Hb ^ s + 'y)/T{Hb — > luX) remains practically unaffected by the non-perturbative 
corrections. 

More interesting effects arise in the photon spectrum: to lowest order it is given 
by a single line reflecting the two-body nature of 6 — *> s + 7; both perturbative and 
non-perturbative corrections turn it into a continuous spectrum. Yet this will not be 
discussed here; the reader is referred to a quickly expanding literature |T3|, ^, Q on 
that topic. 

3.5 Charm Multiplicity in the Final State 

Because of the KM hierarchy |V(c6)p ^ |V(m6)P and since the transition b -^ ccs 
occurs one realizes that 

Number of charm states , , 

Ncharm = ^-^ > 1 41 

B decay 

has to hold where charmonia states enter into the book keeping with a charm multi- 
plicity of two. Since there is not enough phase-space for b -^ cud, civ to transmogrify 
itself into B — > DDD, DDDlv etc., one knows without ado that each of these tran- 
sitions will yield exactly one charm state per B decays. It used to be stated that one 
actually predicts Ncharm = 1-15 because of BR(b — >• ccs) = 0.15. The prediction for 
the latter quantity depends of course quite sensitively on the values adopted for the 
ratio of the charm to the beauty quark mass, as expressed by the function Io{x, x, 0) 
defined in eq. (35a) with x = {mc/nibY. Using values for ttIc and mb as inferred from 
Heavy Quark Expansions one actually finds 

Ncharm ~ 1-2 - 1.3 (42) 

The data exhibit a considerably lower charm content, namely [^ 

A^charm = 0.932 ±0.10 ARGUS (43a) 

Ncharm = 1-026 ± 0.057 CLEO (436) 
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Taking eqs.(42, 43) at face value one would have to state the existence of a significant 
'charm deficit'. On the other hand one should keep in mind that in the transition b — » 
CCS the energy release is not that large. Thus the l/mf, expansion has to be applied 
with a grain of salt here: the non-perturbative as well as perturbative corrections 
in 6 ^ CCS could be larger than expected and negative. But one conclusion can be 
drawn in any event: attempting to lower the predicted semileptonic branching ratio 
for B mesons by increasing BR{b -^ ccs) would fail! 

3.6 Final State Interactions 

Invoking strong final state interactions (FSI) - in the form of phase shifts and ab- 
sorption - to escape conflict with the data represents a time-honoured and popular 
tool in phenomenological analyses; FSI are typically introduced in an ad-hoc fashion 
with an obscure dynamical foundation. Two attitudes towards FSI seem to prevail in 
the community: (a) It constitutes an a priori hopeless enterprise to account for FSI 
even in a semi-quantitative way because they are generated by strong re-scattering 
processes among real hadrons. (b) Whatever the origin of FSI, they are actually 
suppressed in heavy flavour decays. 

Such a situation is quite unsatisfactory, not because FSI are by themselves enlight- 
ening - they certainly are not - but because they represent a conditio sine qua non 
for the observability of direct CP violation in B decays. Attitude (b) would suggest 
that direct CP asymmetries will be small in B decays; while attitude (a) on the other 
hand holds out the hope that one might find sizeable direct CP asymmetries, it would 
make their interpretation quite ambiguous. 

The experience we have gained recently from the QCD treatments of heavy flavour 
decays leads us to the following expectation: FSI in exclusive decays may well pos- 
sess a rather complex and non-trivial pattern that we cannot predict at present; 
yet re-scattering in inclusive transitions must be treatable in QCD using the same 
methodology that has been introduced in the previous sections. 

One estimates FSI effects in non-leptonic B -^ -^charmless decays by analyzing the 
interference of the Penguin amplitude b -^ s[d]qq with the tree- level KM suppressed 
b -^ uus[d] amplitude ||3^. These two amplitudes possess a "hard" FSI phase differ- 



ence 6p due to the on-mass-shell intermediate c and c quarks in the Penguin process; 
one easily finds (5p ~ 0.1. It was shown |^ that all higher order gluon corrections can 



be summed up in a compact way; they reduce 6p by just ~ 20 %. This demonstrates 
that the perturbative expansion makes sense also for this quantity and thus refutes 



earlier claims to the contrary. An elegant analysis can be given ^1| that involves 
rather general arguments based on gauge invariance and the equation of motion in 
QCD, and shows the strong suppression of the decays of the type b ^ s + gg[ggg] ; 
the latter can actually be traced back to the first genuine QCD study of strange 
decays []32| . 



Only the purely perturbative corrections to FSI have been addressed in ref. ||3T 
Non-perturbative contributions still await a detailed analysis; on general grounds one 
expects them to be sizeable, yet at the same time theoretically tractable! 
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3.7 B Decays into Charmonia 

A very intriguing class of reactions is provided by the decays into charmonia states: 
B — >■ ip/ip' /x + X. No rigorous QCD analysis of these modes has been given yet, 
but reasonable phenomenological treatments do exist. Folding the wavefunctions 
for colour singlet cc boundstates with the cc spectrum in 6 — >■ ccs one can make 
predictions about the relative as well as the absolute rates [^ : 

~ 1 : 0.31 : 0.57 : 0.27 (44a) 

Furthermore 

r(6 ^ X2X) = (446) 

holds in this ansatz. These ratios depend of course on the wave functions chosen, but 
are independent of Ci, C2 and Nq- The absolute widths on the other hand depend on 
these parameters 

T{b ^ ^|Jd^rX) ~ (0.42 - 0.45)(C2 + ^c^^Fo (45) 

while being rather insensitive to the values adopted for the quark masses. For the 
part of the ip wave function that is relevant here is probed also in ip ^ e^e~ and can 
therefore be obtained from the data on that electromagnetic decay. In eq.(45) the 
width for directly produced ip was given; including the feed-down from b — > 'ip'/xi +X 
one finds for the total rate for ip production: 

V{b^^ + X)- (0.52 -r 0.56)(c2 + -rrcif^o (46) 

and, therefore, 

5i?(5 ^ 7/- + X) ~ (0.27 ^ 0.4) %. (47a) 

If the 1/Nc term in eq. (46) is omitted one gets, instead, 

5i?(5^^ + X) -(0.83^1.2)%. (476) 

Comparing these predictions with the data from ARGUS and CLEO 

BR{B -> ipX) = (1.09 ± 0.04 ± 0.07)% (48a) 

BR{B -^ ij'X) = (0.30 ± 0.05 ± 0.03)% (486) 

BR{B -^ xiX) = (0.54 ± 0.15 ± 0.14)% (48c) 

one draws three conclusions: 

• The prediction on BR{B —>■ ip'X) vs. BR{B —>■ ipX) is in good agreement with the 
data. 

• For the natural value 1/A^c = 1/3 the improved parton model expression in eq.(46) 
yields a branching ratio that is too low by a factor of ~ 3 — 4. With 1/Nc = the 
observed branching ratio is reproduced, as it is for 1/A^c — 0.45! This means - not 
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surprisingly - that for b -^ ccs transitions the deviations from this model are large 
and not fully implemented through employing a hadronic wavefunction. It remains a 
challenge to theory to identify those corrections. 

• The prediction on BR{B -^ Xi^) vs. BR{B -^ ipX) is ~ 2 a low. This could be 
blamed on a less than optimal choice for the xi wavefunction. As already mentioned 
the relevant part of the ip and ip' wavefunctions is directly calibrated by the data on 
ip, ip' -^ e~^e~; for the Xi wavefunction such a cross check does not (yet) exist and 
one has to rely on a specific ansatz for the interquark potential, which introduces 
additional systematic uncertainties. There exists of course another possibility |37 



the cc pair could be produced in a colour octet configuration which then sheds its 
colour charge through emission of a soft gluon and transforms itself into a colour 
singlet P wave charmonium state. Naive factorization does not hold anymore in the 
presence of this second source of x production in B decays; its weight depends on 
the probability for a colour octet cc to transmogrify itself into a P-wave charmonium 
state which cannot be predicted from first principles at present. On the other hand 
one can extract its size from the requirement to saturate the observed B —>■ xi^ rate; 
subsequently one can then predict the rate for B —>■ X2X. One finds 

V{B ^ xiX) ~ V{B ^ X2X) 

in contrast to eq. (44). 

4 Exclusive Two-Body Decays of Beauty 

A discussion of exclusive non-leptonic decays has to be opened with a note of caution: 
First Theoretical Caveat: The relationship between inclusive and exclusive transi- 
tion rates is nothing short of delicate! 

This piece of common sense can be illustrated by the following example Q] . Con- 
sider the corrections to the decay width of a {Qq) meson that are induced by Weak 
Annihilation, see Fig. 2 

As indicated there are three cuts in the Qq —>■ Qq forward scattering amplitude 
representing different final states, one with an on-shell gluon and the other two with a 
(slightly) off-shell gluon in the form of a qq pair. Summing over these three cuts yields 
an overall correction |Tp that remains finite even in the limit mq/iriQ — >■ [^ p!T[| . 
However a very striking pattern emerges when one considers separately the three 
'exclusive' channels (a), (6) and (c) 0: The contribution (6) constituting the square 
of an amplitude is positive; in the limit mq/niQ — > it is dominated by a term 
+ (mQ/r7ig)|Tp with T denoting a quantity that is regular in the limit mq/mq -^ 0; 
the contributions (a) and (c) on the other hand represent interference terms that, 
taken together, are of the form —{mQ/m?^\T\'^ for mq/mQ -^ so that the sum of 
(a), (6) and (c) possesses a regular limit. We want to draw the following lesson from 
this discussion: a small effect in an overall rate can be due to large cancellations 
among subclasses of decays. The relevance of this statement will become clearer later 
on. 



■''In the real world these three channels can of course not be distinguished; yet this academic model can 
illustrate the relevant point. 
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4.1 Phenomenological Models 

A popular phenomenological model of non-leptonic two-body decays of charm and 



beauty was suggested in ref . |39[ . There are three main ingredients in all models of this 
type: 

(i) One assumes factorization, i.e. one uses (M1M2I J^J^|-D) ~ {Mi\J^\D) ■ (M2I J^|0) 
to describe D -^ M1M2. (ii) One employs one's favourite hadronic wavefunctions 
to compute (Mi|J^|D). Very recently Heavy Quark Symmetry and Chiral Sym- 
metry (for the light quarks) have been incorporated into these wavefunctions |^ . 



(iii) All two-body modes are then expressed in terms of two free fit parameters 
tti and 02 , with ai controlling the 'class I' D^ — > M1M2 and 02 the 'class 
ir D^ -^ M°M2 transitions; both quantities contribute coherently to the 'class III' 
transitions D^ -^ MiM2- The analogous procedure is applied to B decays allowing 
though for al 7^ a^ and 03 7^ 03 • 

With these two free parameters ai2 (and some considerable degree of poetic license 
in invoking strong final state interactions) one obtains a decent fit for the D° and D~^ 
modes (much less so, however, for Dg decays). The situation is rather similar in B 
decays. One has to point out, though, that this success is helped by a considerable 
ambiguity in the estimates of the matrix elements of the type (7r| J^|i?) and by the 
forgiving imprecision in many of the branching ratios measured so far. 

Exclusive decay rates depend sensitively on long-distance dynamics. In the factor- 
ization approximation large distances enter through simple hadronic matrix elements 
{M\Jfj_\M) and (M| J^|0). The coefficients ai and 02 are: 

ai = ci + ^C2 (49a) 

a2 = C2 + ^ci (496) 

where Ci and C2, see eq. (17), include the radiative QCD corrections due to hard 
gluons; the quantity ^ is introduced as a fudge factor. Literally speaking, factorization 
implies ^ = 1/A^c = 1/3- Quite often ^ is treated as a fit parameter refiecting our 
ignorance of how a quark and an anti-quark that are not correlated in colour combine 
to form a meson. Thus, deviations of C, from 1/3 parametrize, in a certain way, 
non-factorizable contributions. The non-factorizable contributions show up even in 
perturbation theory, but these seem to be numerically small. An example of the non- 
perturbative contribution which can violate factorization is provided by the strong 
final state interactions (FSI). The latter can easily change the value of C,- Some 
prominent features of FSI are actually added in an ad-hoc fashion; yet even so one 
would a priori not expect that all these non-perturbative corrections can by and large 
be lumped into a single fudge factor C,- At this point we would like to express another 
note of caution: 

Second Theoretical Caveat: There exists no general proof of the Dogma of Factor- 
ization for real hadrons; we actually consider it unlikely to he of universal validity, 
in particular for class II and III transitions. It thus makes eminent sense to subject 
this dogma to as many different experimental tests as possible. 

A fit to eleven exclusive hadronic B branching ratios yields ^ 



ai = 1.05 ±0.03 ±0.10 (50a) 
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as = 0.227 ±0.012 ±0.022 (506) 

A few comments are in order here: 

(i) The results stated in eqs. (50) are bad news for a popular program to infer ai 
and a2 from perturbative QCD plus factorization: the effective transition operators 
are renormalized with coefficients c_|_,c_ generated from perturbative QCD; then the 
hadronic matrix elements are factorized. The coefficient 02 does not come out cor- 
rectly. Moreover, the so-called rule of 1/Nc (this is an ad hoc assumption that only 
the terms leading in 1/Nc are to be retained), which miraculously helped in the two- 
body decays of charm, seems only to worsen the situation for beauty. The program 
"factorization plus the rule of 1/Nc" implies ai = ci, 02 = C2, in clear conflict with 
the findings in B decays, eqs. (50), since C2 is negative! Yet even before the following 
statement applied: 

Third Theoretical Caveat: It is very unlikely that the rule of retaining only leading 
terms in 1/Nc is universally implemented in QCD. Our analysis of inclusive heavy- 
flavour decays actually found cases where this rule was dynamically realized (a) for 
D as well as for B decays, (b) for D, but not for B decays, or (c) for neither [|, ^, ^ 
(see also Sect. 4.2). 

(ii) The most striking feature of eqs. (50) is that the relative sign between ai 
and 02 is positive, i.e. that a constructive interference occurs in class III transitions, 
i.e. in exclusive two-body B^ decays! This is surprising in four aspects: (a) It is 
in clear contrast to the situation in D decays where the D~^ modes suffer from a 
destructive interference, (b) While it does not pose any fundamental problem for 
the BSW model, it represents a basic failing for the 1/Nc ansatz which predicts a 
destructive interference in the two-body modes of B as well as of D decays, (c) The 



observed enhancement of the B~^ rates is remarkably uniform pl| : 

BR{B- -^ D\-)/BR{B'^ -^ D+tt") = 1.71 ± 0.38 (51a) 

BR{B- -^ D^p-)/BR{B^ -^ D+ p-) = 1.60 ± 0.46 (516) 

BR{B- -^ D*^TT-)/BR{B'^ -^ D*+7r") = 1.79 ± 0.39 (51c) 

(d) It raises the question of whether the same constructive interference might occur 
for the inclusive rate thus shortening t{B~) relative to T{Bii) rather than lengthening 
it. 

As stated in Sect. 3 the B~ lifetime is predicted to exceed the Bd lifetime, albeit 
by a few percent only; yet one has to keep in mind the First Theoretical Caveat stated 
in the beginning of this section: a small correction in the inclusive rate is quite likely 
to be made up by large contributions of alternating signs coming from different classes 
of exclusive transitions. 

For a better understanding of these problems one has to progress to treatments 
that are rooted more firmly in QCD. However we would like to first stress the im- 
portant lessons we have learnt and are still learning from these phenomenological 
descriptions: 

- They have yielded quite a few successful 'predictions' in a user-friendly way 
providing a useful pattern for cataloguing dozens of the decay rates, a reference frame 
for more sophisticated approaches. 
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- They have helped us considerably in focusing on the underlying theoretical 
problems such as the question of factorization or the 1/Nc rule. 

- From their fits to the data they provide us with valuable, albeit indirect infor- 
mation on the final state interactions. Such information is crucial in studies of direct 
CP violation. 

4.2 Theoretical Analysis 

The first treatment of two-body decays of heavy flavour hadrons that is intrinsically 
connected to QCD was given in ref. |Q some time ago for D decays, within the 



framework of the QCD sum rules. To extend this analysis to non-leptonic B decays 
is a rather non-trivial undertaking; work on treating B -^ J /ipK through QCD is in 
progress now p4 



Recently there has been progress in analysing deviations from the 1/Nc rule 
(see Sect. 4.1) in some exclusive two-particle B decays |^. It has been shown 
how non-perturbative effects in QCD can provide a dynamical realization of this 
rule in some decay channels, but not in others. The key role is played by the same 
chromomagnetic operator a ■ G that was repeatedly discussed above in connection 
with the inclusive decays. It generates non-perturbative corrections that are specific 
for the decay channel considered, yet can be estimated in a model-independant way. 

The general method can be illustrated through the example of the decay B^ —>■ 
D^Ti^. The non-factorizable terms in this decay are reducible to the amplitude 
(DnlCoctlB) where 

Coct ~ (c7m(1 - 75)t"6)(rf"7M(l - 75)t"M), 

with t"" denoting the colour matrices. Colour has to be exchanged between the brack- 
ets because otherwise the light quarks can not form the pion; this can be done by a 
soft gluon. 

To estimate this effect one considers the time ordered correlation function 

^/3 ^ j d^xe''^^D\T{Coct{x),Af}\B) (52) 

where an auxiliary axial current A^ = u'^^'^r^d annihilates the pion and q is an external 
momentum flowing through A^ . To calculate the correlator one adopts a similar 
procedure as in QCD sum rules. After continuing A^ into the Euclidean region 
—q^ = Q^ ~ 1 GeV ^ one invokes duality in the following way: on the one hand one 
expresses A^ in terms of Mnon-fact, the non-factorizable part of the amplitude, as 
obtained from Cod'- 

^°'=*(g2) = Mr,on-fact^ + -. (53) 

where -|-... denotes the contributions from the higher resonances produced by the axial 
current A^; on the other hand one applies an OPE to A^ to find 

^°'\Q') = ^^,^{D\cJ,l5Go.,b\B) + ... (54) 
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where +... now denotes pre-asymptotic corrections from higher dimensional operators. 
If both masses, rric and rub, are treated as heavy, the matrix element in eq. (54) is 
related to the expectation value of the chromomagnetic operator aG and is expressed 
in terms of (Iq, see eq. (36b). Comparing the leading terms in eqs. (53, 54) one 
obtains for the ratio of the non-factorizable to the 1/Nc factorizable parts of the 
amplitude 

'non-fact . 2/2 V / 

An essential difference to the usual QCD sum rules analysis is worth noting: rather 
than the vacuum condensates we deal here with hadronic expectation values between 
heavy flavour hadron states. Numerically one finds 

rnon-factiB"^ ^ D+Tl-) ~ -1 (56) 

and likewise for the mode B^ -^ D^p^; i.e. in these two modes non-factorizable con- 
tributions which are necessarily of order I/Nq basically cancel against the non-leading 
factorizable contributions of order 1/Nc thus leading to a dynamical realization of the 
1/Nc rule. The findings are similar for B -^ DD decays. Yet the situation changes 
for other modes: in B ^ D*D* the non-factorizable contributions are suppressed; 
likewise in B° -^ D*+7r-, D*+p~. 

There are two general conclusions we want to draw from this analysis: (i) The 
weight of non-factorizable contributions is indeed quite channel dependent - as ex- 
pected, (ii) The discussion so far was given for class I transitions where we found 
that the coefficient ai is not quite universal, but has some channel dependence. The 
situation is much more complex and actually at present unclear for class II (and III) 
transitions. For the 02 amplitude contains the matrix element {B\b'~f^gGauu\TT) ■ This 
formfactor, unfortunately, can not be determined from HQET. 

This treatment has reached so far only the qualitative or at best semi-quantitative 
stage. No detailed analysis has yet been given about the question whether one can 
really suppress the contributions from higher resonances and from higher dimensional 
operators in the sum rule of eqs. (53, 54) to a sufficient degree. Furthermore, the 
operator product expansion, eq. (54), assumes that {Mb — Md)/{Mb + Mjj) « 1; 
in the real world this parameter is rather of order unity. Also radiative corrections 



have not been included yet. Further details can be found in ref. |^^. Thus, much 



more theoretical work is needed before such a treatment finds its definite form. 

4.3 Prizes to be Attained 

The theoretical methods one applies to exclusive decays are often not of the most 
lucid kind. Yet they are essential (if imperfect) tools for addressing fundamental 
questions. Let us cite just one topical example: 

AT{Bs), i.e. the lifetime difference between Bg^short and Bg^iong, is usually com- 
puted from the quark box diagram with internal c (and u) quarks, leading to a result 
like the one quoted in eq. (15). However the weight of such a short- distance contribu- 
tion to Ar{Bs) is much more uncertain than that of the local contribution to Am{Bs), 
which is given by virtual top exchanges. It is therefore conceivable in principle that 

21 



nontrivial long-distance dynamics could provide in reality a larger contribution. This 
possibility has been analysed by the Orsay group |H6[ considering separate transitions 



and it has been inferred from present data Ar/T{Bs) ~ 0.15 in satisfactory agreement 
with what one finds using the quark box diagram, eq. (24). 

5 Summary and Outlook 

5.1 Status 

New and more powerful second-generation theoretical technologies are emerging: 
QCD sum rules, Heavy Quark Symmetry, l/mq expansions and lattice simulations 
of QCD. They are leading to 

• significant conceptual progress, namely a better understanding of (i) the form 
and size of non-perturbative corrections, (ii) the relationship between charm and 
beauty decays, where the former play the role of a microscope for the non-perturbative 
corrections in the latter, and (iii) the differences and similarities of baryon vs. meson 
decays; 

• the realization that charm and beauty baryons deserve detailed studies in their 
own right. 

• a quantitative phenomenology that is genuinely based on QCD: 

^^^ ^ 2; BRsl{D+) r-. 16%, 5i?5Lp°) ~ 8%; 






1.0 ± a few per cent , (57) 



i.e. the data are reproduced within the accuracy of the expansion. 

7ll^-^ + °-°^- (200MeV)^ ^ BRs.iB)> 12% (58) 

^£(^^018. f-^ I^~09 (59) 

V{B,) -^-'^ (200MeV)2' r{B,) ^'^ ^'^^ 

dV{B- -^ lvX)/dEi ^ dV{Bd -^ lvX)/dEi (60) 

The prediction for BRsl{B) is somewhat larger than present CLEO and ARGUS 
measurements. This could turn out to be a serious - or intriguing - discrepancy. It 
could conceivably signal the presence of anomalously large higher-order contributions 
that so far have not been included in the theoretical expression. In that case one 
would expect lifetime ratios for beauty hadrons to differ more from unity than stated 
in eqs. (58) and (59) p^. 

Phenomenological models for nonleptonic two-body modes are encountering dis- 
crepancies with more precise data; yet they continue to be useful and help us in 
focusing on the underlying theoretical issues. 
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5.2 Future 

One can expect a refinement of and increased cooperation (rather than just coexis- 
tence) between the second-generation theoretical technologies. On the experimental 
side one can hope for 

- lifetime measurements for individual beauty hadrons with 10% accuracy soon 
and percent accuracy in the longer run; 

- data on t{Bs) separately from Bg -^ ipcj) and from Eg -^ IvOg] 

- perform the 'class I, II, III' phenomenology individually for iTM allowed and KM 
suppressed B decays. 

The primary goal in all these efforts is to be able to exploit to the fullest over the 
next 20 years or so, the discovery potential or even discovery guarantee that awaits 
us in beauty physics. It certainly would be a crime not to make all conceivable efforts 
to obtain the required experimental facilities. 
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FIGURE CAPTIONS 

Fig. 1 

Corrections to tlie bq -^ bq forward scattering amplitude induced by Weak Annihila- 
tion. The three cuts (a), (b) and (c) represent different final states for the bq decay. 

Fig. 2 

(a) Diagram generating the operator bb. 

(b) Diagram describing Weak Annihilation. 

(c) Diagram describing Pauli Interference. 
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